Non-Boltzmann classical correction to the velocity auto-correlation function for 

isotropic scattering in two dimensions 
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The classical correction to the velocity auto-correlation function of non-interacting particles due 
to memory effects, which are beyond the Boltzmann equation, is calculated both analytically and 
numerically for the case of isotropic scattering in two dimensions. 
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1. Introduction 

Some fourty years ago a number of pioneering theoretical works appeared [1-7], which were devoted to the calculation 
of the memory-effect (non-Markovian) corrections to the kinetic coefficients and to the velocity auto-correlation 
function, through which they can be expressed. It was shown that such quantities as the diffusion coefficient, the 
electrical conductivity, viscosity etc, can not be expanded in powers of density, as it was assumed previously. Later, it 
was also demonstrated that the velocity auto-correlation function (which we abreviate below as " correlation function" ) 
does not decay exponentially at large times, as the Boltzmann equation predicts, but rather contains a slow power-law 
tail [8-10]. 

These important deviations from the conventional picture based on the Boltzmann equation are due to memory 
effects neglected in the Boltzmann approach. For the case of non-interacting particles performing diffusion in a 
static random potential field these effects are due to returns of the particle to previously visited regions. In the 
two-dimensional case it was shown by Ernst and Wejland [9] that, because of returns, for times t much greater than 
the mean free flight time r, the correlation function has a negative tail decaying as 1/t 2 . Bruin [10] has performed 
numerical calculations for scattering by hard disks, which showed that this asymptotic behavior for t >> r may 
appear only at very long times, it is not reached even at t ~ 10r. 

Returns after a single collision (see Fig. la) are important at t < r. For t << t they give a contribution to the 
probability of return increasing as 1/t. This, in turn, leads to a non-analytical in the small parameter d/£ = Nd 2 
correction to the diffusion coefficient on the order of (d/£)ln(£/d), where d is the effective scattering diameter, 
I = (Nd) -1 is the mean free path, and N is the two-dimensional concentration of scatterers [2,5,11,12]. Returns after 
two or more collisions (see Fig. lb) give a smaller correction on the order of d/£. 





X 



FIG. 1. Illustration of the memory effect due to return to the same scattering center after a single collision (a) or several 
collisions (6) 



In our previous work [13] we have drawn attention to another memory effect (the corridor effect) important in 
backscattering events. If a particle travels a distance x after which it is backscaterred and returns to the initial point, 
the probability of this round trip of length 2x is proportional to exp(— x/l), not to exp(— 2x/£), as would suggest 
the conventional Boltzmann approach, since the existence of a free corridor of width d allowing the first part of the 
journey garantees a collisionless return. This effect also gives a correction to the diffusion coefficient on the order of 
d/L 
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The role of memory affects becomes greatly enhanced in the presence of an applied magnetic field since the circling 
motion of an electron naturally increases the probability of returns, and many works were devoted to classical magne- 
toresistance, which for a degenerate 2D electron gas is entirely due to memory effects [14]. In particular, it was shown 
that the irreversibility introduced by the magnetic field destroys the corridor effect resulting in an anomalous magne- 
toresistance in classically weak magnetic fields [13]. An analytical theory of the anomalous low-field magnetoresistance 
was recently developped in Ref. 15. 

The existing derivation [6,7] of the memory effect corrections is rather long and cumbersome being based on a 
special technique of the so-called ring diagrams. 

In the present paper we develop a relatively simple approach to the problem. We consider classical non-interacting 
particles with a fixed energy in two-dimensions scattered by randomly positioned centers with a given differential 
cross-section in the absence of magnetic field. We start with deriving a modified kinetic equation, which takes into 
account single returns after an arbitrary number of collisions. We then derive the leading logarithmic correction to the 
diffusion coefficient for an arbitrary angular dependence of the differential cross-section. Finally we find analytically 
the complete time dependence of the velocity correlation function for the special case of isotropic scattering taking 
into account the corridor effect and calculate the correction to the diffusion coefficient for this case. We perform 
numerical simulations for the cases of isotropic scattering and scattering by hard disks, and we find a good agreement 
between the numerical and analytical results. 

2. Derivation of the basic equation 

The derivation of the conventional Boltzmann equation from the Liouvillc equation for a particle moving in the 
presence of randomly positioned scattering centers involves a number of simplifying assumptions. First, it is assumed 
that the mean free path, £, is much greater than the effective scattering diameter d or, otherwise that Nd 2 << 1, 
where N is the concentration of scatterers. Second, it is assumed that the distribution function does not change much 
on the space scale on the order of d and on the time scale on the order of d/v, where v is the modulus of the particle 
velocity, which is conserved since the scattering is elastic. These two conditions make it possible to describe scattering 
in terms of the differential cross-section a{4>) attributed to scattering centers positioned at given points in space. 
Third, memory effects are neglected: the Boltzmann approach is equivalent to randomly redistributing the scattering 
centers after each collision. Neglecting memory effects allows to take the average over the positions of scatterers in 
the Liouville equation, rather than use the average of its solution. 

Restricting ourselves by the first two assumptions only, we can write down the following equation for the distribution 
function /(r, </>, t) (see Appendix A and Ref. 15 for the derivation of this equation in the simplest case of scattering 
by hard disks). 



Note, that the distribution function, /, in Eq. (1) depends on the positions of scatterers, rj, for a given realization. 
Also, since the scattering is elastic, the distribution function depends only on the polar angle, 0, of the velocity vector 
v. The conventional Boltzmann equation can be obtained from Eq. (1) by replacing the actual density of scatterers, 
J2i <K r ~ r i), by its average value, N. 

The correlation function, K(t), (normalized by the condition K(0) = 1) and the diffusion coefficient, D, can be 
expressed through the solution of Eq. (1) with the following initial and normalization conditions: 




(1) 



where Yi are the coordinates of the scattering centers, and t is the scattering operator proportional to the Boltzmann 
collision integral: 




(2) 






(3) 



Then 
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(4) 




(5) 
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The angular brackets in Eq. (4) denote averaging over the positions of the scattering centers. 

We now develop a method to calculate K(t) and D taking into account non-Boltzmann memory effects due to 
returns of the particle to previously visited scattering centers. As a first approximation we choose the solution 
G(r — r',(j>,(j)' , t) of the conventional Boltzmann equation: 

dG dG A ~ „ 

- +V - + vNTG = (6) 

with the initial condition 

G{v-v',ci>, < p',Q) = 5{v-v l )6{<t>-<t>'). 

We re- write Eq. (1) in the form 

%+v%+Nvff = -v(r)ff, (7) 

where ^(r) = J^. <5(r — rj) — N is the fluctuation of the scatterer's concentration, and (f(r)) = 0. As above, the 
distrubution function / in this equation depends on the actual positions of the scattering centers rj. 
Eq. (7) may be also written in an integral form using the formal solution of Eq. (6): 

f(r,<f>,t)=G(r,<f>,0,t)-vGvff, (8) 

where G is the integral operator with the kernel G(r — r', cf>, <fi' , t — t'). 

We now substitute / given by Eq. (8) into the right-hand side of Eq. (7) and we take the average of the resulting 
equation over the positions of the scatterers rj. In doing this we must deal with the product of functions is(r), ^(r'), 
and the distribution function /, which all depend on the coordinates r 4 . We decouple / from the averaging procedure 
by writing 

MrMr')/> « <i/(rMr')></>. 

This approximation takes into account single returns to the same scattering center, but neglects multiple returns. 
In the absence of correlation in the positions of the scatterers, which we assume to be true, we have 

{v{r)v{r')) = N5{v-t'). 

Hence, we obtain the following equation for the averaged distribution function (to simplify the notations we replace 
hereafter (/) by /): 



dt 



+ w^J- + Nvff = Nv 2 [ dtp [ dt'TG{0, <j>" -<f/,t- f)ff. 
dr Jo Jo 



(9) 



In the right-hand side of this equation the function Tf is a function of 0', the left operator T acts on a function of 
the variable 4>" . 

The Green function G(0, <f>" — 4>',t) is the probability for a particle to return to the initial point after time t with 
velocity directed at angle <f>", provided that the initial velocity is directed at angle <f>'. Obviously, this probability 
should depend on the difference <p" — <p' only (this is not the case for G(r — r', ((>", </>', t) if r ^ r' ). 

The correction in the right-hand side of Eq. (9) coincides with the result of Weijland and Van Leeuwen [7] obtained 
by the ring-diagram technique. 

Eq. (9) is the basis for the following calculations. The solution of this equation is needed for calculating the 
correlation function K(t) and the diffusion coefficient D with the help of Eqs. (4), (5). In fact, because of the 
approximations made during our derivation the right-hand side of this equation should be considered as a small 
perturbation. 

3. Derivation of the correlation function and the diffusion coefficient 

To calculate K(t) we need to know the integral of the distribution function over the coordinate: 

F{cj>,t) - J f(r,<j>,t)dv. (10) 
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Integrating Eq. (9) over r we get 



dF 



2tt ft 



^- + NvTF = Nv 2 l d<t>' I dt'TG{Q,(j)" -<f>', t-t')TF. (11) 



This equation can be further simplified by using the property that the functions exp(im,(j)) are the eigenfunctions 
of the scattering operator T: 

NvT exp(im(f)) — 7 m exp(im^), 

7m = Nv / (1 - cos(m0))cr(0)#. (12) 
Jo 

Expanding the functions F and G in Eq. (11) in Fourrier series, we obtain 

' IF '" l m F m = ^r7lJ*dt'G m (0,t-t')F m (t'). (13) 



dt 
where 



1 f 2lT 

Fm ^ = 2^J c ^v{-ifn(j))F((j),t)d(j), 
1 f 2n 

Gm(0,t) = — exp(-im<p)G(0,4>,t)d(t>. 
27T Jo 



Eq. (13) allows to calculate any moment of the distribution function with the corrections due to returns taken into 
account. 

As seen from Eqs. (4), (5), and (10), the correlation function and the diffusion coefficients can be expressed via the 
functions Fi(t) and -F-i(i) = Fi(t): 

K(t) = 2irReF 1 (t). (14) 

Eq. (13) can be solved exactly. However, it should be noted that this equation, like our basic Eq. (9) takes 
into account single returns only and consequently, it gives the correct results only in the leading order in the small 
parameter Nd 2 . For this reason, as we have already mentioned above, the right-hand term in Eq. (13) should be 
considered as a small perturbation. This means that in this term one can use the Boltzmann expression for F m (t): 

F m (t) = exp(-7 m i), (15) 

where the factor (l/27r) appears because of the initial condition F(</>, 0) = S((f>). Substituting this expression onto the 
right-hand side of Eq. (13), for m — 1 we get: 



dt 
where 



' lFl ^~nIo Gl(0 ' 1 ~ t] exp (~ 7 *'^'' ( 16 ) 



,-2-k 

7 = 71 = t _1 = Nva tr = Nv / f(0)(l — cos(f))d(j) 

Jo 



is the inverse momentum relaxation time, or transport time, r. 

The solution of Eq. (16) can be easily found. After some manipulation it can be presented as follows: 



F 1 (t) = -L exp(-Tt) + ^ J Gi (0,t - f) cM-lt')t'dt'. (17) 
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Finally, using Eq. (14) we obtain the following expression for the correlation function: 

K(t) = K (t) + SK(t), SK(t) = 2tt^- / Gi(0, t - t') exp(-jt')t'dt' , (18) 

N Jo 

where 5K(t) is the correction to the Botzmann result, Ko(t) = exp(— "ft). 

Substituting this expression into Eq. (5) and changing the order of integrations over t and t', we obtain a formula 
for the diffusion coefficient: 

2 2 poo 

D = D + SD, A> = ^, SD = ^ r j^ Gi(0,t)dt. (19) 

Here Dq represents the Boltzmann result, while 6D gives the correction due to returns, which is beyond the 
Boltzmann equation. One can see that calculation of these corrections is reduced to finding G\ (0, t) , which is equivalent 
to finding G(0, <f), t), the Boltzmann probability of return to the initial point r = with velocity directed at angle 4> 
with respect to the initial velocity. 



4. Contribution of returns after a single collision 

We have mentioned in the Introduction that the leading (logarithmic) correction to the diffusion coefficient is due to 
returns after a single collision (see Fig. 1). We will now show that the contribution of this process can be easily found 
for an arbitrary angular dependence of the scattering cross-section a{4>). To do this, we must find the corresponding 
contribution to the function G(0, </>, t). 

We re-write Eq. (6) in an integral form: 

G(r, <f>, t) = S(r - v t)5(<t>) cxp(- la t) 



+Nv d(/)'<j((f)-(f)') exp(-j (t-t'))G(r-v'{t-t'),(j)',t')dt', (20) 
Jo Jo 

where 70 = Nva, a is the total scattering cross-section, vo is the initial velocity (with <j> = 0), and <j>' is the angle of 
the vector v'. 

Eq. (20) is convinient for obtaining the contribution of a given number of collisions before return by iterations. 
Since we are now interested in returns after a single collision, we can insert into the right-hand side of Eq, (20) the 
zero approximation for G, equal to S(r — v t)5(<j)) exp(—j a t). Denoting the correction thus obtained as SG, we get: 

P 27T f>t 

5G(r,<f>,t) = Nv dcf>'a((f>-(j>') exp[- 7o (i - if) - -yot']S(r - v'(t - t') - v t')dt'. (21) 
Jo Jo 

Putting r = and taking into account the relation 

f)+vo*0= Wz£j*( t -2f), 



we find 



Thus, for SGi(0,t) we obtain: 



5G{0, <t>,t) = — 6XP( 5(<f> - 7r)(r(0). (22) 



To find the correction SD to the diffusion coefficient we must substitute this expression in Eq. (19). The integral 
in Eq. (19) diverges logarithmically, which is the result of ignoring the finite radius of the scattering center. To avoid 
the divergency, we will replace t in the denominator of the expression Eq. (23) by t + t 0} where the cut-off time to 
is on the order of a /v. Such a regularization seems reasonable, since it makes the return probability to be finite at 
t — > 0, as it should be if the finite radius of the scattering center is taken into account, and since our entire theory is 
applicable for time scales larger than a /v. 
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Eq. (23) must be also corrected to include the corridor effect [13,15], which we mentioned in the introduction and 
which is not taken into account by Eq. (9). This effect is relevant for backscattcring events, when the particle follows 
practically the same path (in the opposite direction) after a collision. The probability to make this round trip without 
collisions during time t should be exp(— jot/ 2) rather than exp(— jot) as Eq. (23) says, because once the path 1 — ► 2 
during the time t/2 is collisionless, we are sure to have no collisions on the return path 2 — > 1. This effect can be 
taken into account by a more accurate evaluation of Eq. (1), which (within the approximation of point-like scatterers) 
contains all memory effects, including the corridor effect, see Ref. 15 where this was done for scattering by hard disks. 
Here we will simply modify Eq. (23) "by hand" replacing exp(— j t) by exp(— jot/2). 

Thus, we replace Eq. (23) by 

6G 1 (0,t) = -^ eM -^ t/2) a(n). (23a) 
v ; 2ttv t + t 

We then obtain the correction to the diffusion coefficient due to returns after a single collision: 

S -^ = Na tr a(TT)(-\n(^-) + C-\n2), (24) 
L>o \ 70*0 1 

where Dq = v 2 /2j is the Boltzmann value of the diffusion coefficient and C = 0.577 is the Euler constant. Since 
the exact numerical coefficient in the argument of the logarithm is unknown, the constant C — In 2 could be safely 
discarded. We prefer to keep it in order to compare Eq. (24) with the more general formula derived in Sect. 5 and 
to have a clear definition of the cutoff parameter to when discussing numerical results in Sect. 6. 

With regard to this formula it should be noted that the argument of the logarithm contains the constant 70, which is 
expressed through the total scattering cross-section. The latter is equal to infinity for any realistic scattering potential 
that does not drop to zero at a finite distance (as it is the case for scattering on hard disks when the total scattering 
diameter is of course equal to the disk diameter). The divergency of the total cross-section is due to very small 
scattering angles. Since we consider the scattering centers as points, any scattering event will deviate the particle 
from its path connecting two scattering centers, however small the scattering angle may be. If the finite radius of the 
scatterer is taken into account, then scattering angles less than ~ d/£ will not matter anymore. This means that the 
integral over <j), which gives the total scattering cross-section, entering the definition of 70, should in fact be truncated 
to exclude scattering angles less that d/l. Because of the logarithmic dependence of SD on 70, the exact value, which 
should be attributed to this cut-off, is not very important. 

Substituting Eq. (23a) into Eq. (18), we obtain the correction to the correlation function in the time interval 
t « t « 7 f 7 1 : 

6K(t) = -Na tr a(ir)jtln(t/to). (25) 

To conclude this section we note that returns after multiple collisions, which are important for times t > r, give 
a correction to the diffusion coefficient on the order of Na 2 r . Thus, as one can see from Eq. (24) the logarithmic 
correction dominates if a(n) is not too small, i.e. in all cases when the backscattering is not strongly supressed. 



5. Correlation function for isotropic scattering 

For the special case of isotropic scattering a tr — &, cr(n) — &/2ir, and 70 = 7 = r _1 = Nva. This simplification 
allows to find analytically the Green function G(0, <j>, t) for arbitrary times (that are larger than the cutoff time to). 

The details of the calculation are given in Appendix B. The function Gi(0, t), through which the corrections to the 
correlation function and the diffusion coefficient are expressed, has the form: 

Gl (0, t) = (- + 1-°*P(-*) + 2 exp(-x) (Ei(x) - Ei(2x)) ) , (26) 

where we have introduced the notations x — t/r, x = t /r << 1, and Ei(a;) = — J^° x dyexp(—y)/y is the integral 
exponent. Since to ~ &/v an d t = (-/Vctv)" 1 , the cutoff parameter can be presented as xq — aNa 2 , where a is an 
unknown constant on the order of unity. 

Comparing Eqs. (26) and (23a), one can see that the first term in Eq. (26) represents the contribution from returns 
after a single collision with the cutoff at t — > and the corridor effect taken into account as described in the previous 
section. The other terms in Eq. (26) correspond to returns after two or more collisions and do not contain any 
singularities. 
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A direct calculation using Eqs. (19) and (26) gives us the correction to the diffusion coefficient: 



D 2tt V V V > 



(27) 



which differs from Eq. (24) describing the contribution of returns after a single collision by an additional term (— In 2) 
in the brackets, which is the contribution of returns after multiple collisions. 

The correction to the correlation function, 8K{t) can now be calculated by using the Eqs. (26) and (18). A direct 
calculation gives the following, rather cumbersome, expression (here again x = t/r): 

Act 2 

SK(t) = 

k(x) = 3e- x - 2e- x ' 2 - 1 + 2e~ 2x (Ei(2x) - Ei(x)) + (x + x )e~ x (Ei(^t^) - Ei(y)) 



+e~ x 



(a;-2)(Ei(a;) - C-lnz) + x + 2(x - 1) ln2 



(28) 




10" 
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10 10 10 
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FIG. 2. a) Plot of k(x) — —(2ir/Na 2 )SK(x) versus x = t/r for different values of the cut-off parameter xq. Continuous, 
dotted and dashed curves correspond to x = 0.01, 0.02 and 0.04, respectively. The inset shows the dependance of the maximum 
value, km, on ma^o- b) The long time part is emphasized in a log- log plot to show that the different curves become quickly 
superimposed for x > 5 and reach very slowly the asymptotic regime, l/2a; 2 , as x — > co (shown by the long-dashed line). 

The function k(x) is presented in Fig. 2a for several values of the cutoff parameter xo- In inset is shown the linear 
dependance of the maximum value of k(x) on \nx n . 
For small times, xq « x « 1, Eq. (28) gives 

k(x) = xln(a;/a;o), (29) 
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so that SK(t) coincides with the expression given by Eq. (25), if one takes into account that for isotropic scattering 
&tr — & and (t(tt) = a/2n. 

At large times, x >> 1, the correlation function is dominated by the contribution of returns over long diffusive 
trajectories involving many collisions, ft can be presented as a series in inverse powers of t. Using Eq. (28) in this 
limit, we obtain the leading terms of this expansion: 



(30) 



The leading term (<~ 1/t 2 ) coincides with the known result obtained by Ernst and Weyland, [9] while the second 
one gives the first correction. We note that the leading term (but not the correction) is a universal result, which can 
be derived in the diffusion approximation [9] and which does not depend on the scattering cross-section. It can be 
seen that the asymptotic 1/t 2 behavior is approached very slowly, so that even for t = 50r (x = 50) the correction 
still makes 10% (see Fig. 2b). This explains why the predicted asymptotic behavior was not reached in the numerical 
simulations done by Bruin [10], as well as in our simulations to be presented in the next section. 




x=t/x 

FIG. 3. Plot of k{x) (continuous line) and k(x) (dot-dashed line), given by Eqs. (28) and (B8), respectively, for the same 
value of the cut-off parameter, xq = 0.01. The relative difference (k — R)/k is shown in the inset. 



Finally, in Fig. 3, we present, for comparison, the function k(x) given in Appendix B, Eq. (B8), in which the 
corridor effect is not taken into account. One can see that the role of the corridor effect is noticable, but rather small. 



6. Numerical simulations 

Using a random number generator, the centers of AT disks of diameter d are uniformly randomly positionncd on a 
plane inside a square box of edge length L. We take L/d = 1000 to be sure that L remains more than one order of 
magnitude larger than the mean free path, even for the smallest concentrations that we have considered. The value of 
Af is chosen to obtain the desired value for the dimensionless concentration, Afd 2 /L 2 = Nd 2 . Numerical simulations 
were performed for Nd 2 equal to 0.064, 0.032 and 0.016. To calculate the velocity correlation function K(t), as well 
as the diffusion coefficient D given by Eq. (5), we first choose an initial point at random inside the box with an 
initial velocity direction arbitrarily chosen along the x-direction. We then determine the trajectory of a point-like 
particle by joining the successive impact points, an impact point being the first intersection, calculated analytically, 
of the linear trajectory with a disk periphery, the particle always coming from outside of the disk. We use standard 
numerical tricks to accelerate the search for impacts. When choosing the initial position, the disk interiors are not 
excluded (it was checked that excluding them introduces only a weak numerical difference, which vanishes in the 
limit Nd 2 — ► 0). The trajectory is made of successive straightline segments between collisions. Periodic boundary 
conditions are imposed at the edges of the square box. 

In the hard-disk (or Lorentz) model the scattering angle <f) is related to the impact parameter p by p — (d/2) cos(0/2) 
(we consider < 4> < 2tt, so that p may be negative). The differential cross-section <r(<p) — \dp/d(j>\ = (d/4) sin(0/2) 
is anisotropic, with enhanced backscattering (in contrast to isotropic scattering by hard spheres in three-dimensions). 
As a consequence, the transport cross-section is larger than the total cross-section: a tr — (4/3)d. 

To simulate isotropic scattering we define the scattering angle by the relation p = (d/2) (1 — (fr/ir), so that o{4>) = 
d/2ir. This relation is somewhat artificial, it leads to unrealistic situations, in which the disk may be cut by the 
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scattering trajectory. This isotropic model, which is built to have an angle-independent differential scattering cross- 
section, does not correspond to any realistic potential and is considered for the sole purpose to check our analytical 
theory. 

In order to calculate the correlation function, K{x), and its dimcnsionless integral up to a time t = xt, I(t) = 
J^ T K(x)dx = X(t)/vr, where X(t) is the overal displacement of the particle in the x-direction during time i, we 
consider a set of 5000 discrete values x n of the reduced time x = t/r, regularly spaced up to x m — 20, and we calculate 
both quantities for each x n . The values of K{x n ) and I(x n ) are averaged over 100 disk configurations and 10 s trials 
for the starting point, except for the lowest concentration Nd 2 =0.Q16 where only 10 7 trials for the starting point 
were considered. In this way we get practically continuous curves for both K{t) and I{t). The reduced value of the 
diffusion coefficient, D/Dq, is obtained by extrapolation of the results for I{t) to t — > oo. 




-0.2 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 ■ 1 

2 4 6 8 10 12 

X=t/X 




-1 - 



4 8 12 16 20 

X=t/X 

FIG. 4. a) Results of the numerical calculations for the isotropic model as a plot of K(x) for different concentrations Dotted, 
dashed and long-dashed curves correspond to iVd 2 =0.064, 0.032 and 0.016, respectively. The pure exponential Ko = exp(— x) 
is shown by the continuous curve, b) For t > 4r the scale is enlarged by a factor 10 4 to better show the difference between 
K(x) and Ko(x). 

In Fig. 4, we present the numerical results for K(t) for the isotropic scattering together with the pure exponential 
curve K (t) = e~*/ T expected from the Boltzmann approach. The enlarged scale in Fig. 4b allows to clearly see the 
negative departure from the exponential behavior and the long-time tail of the correlation function. We have also 
calculated numerically the correlation function for the hard disk scattering. The results are qualitatively similar to 
the isotropic case and agree with the numerical simulations of Bruin [10]. 

Figure 5 presents the comparison of the numerical results for the normalized correction to the correlation function, 
k(x), for isotropic scattering with the theoretical formula given by Eq. (28). The only fitting parameter is the the 
constant a in the definition of the cutoff value of xo — aNd 2 . By choosing a — 0.424 we obtain a practically perfect 
fit to the theoretical curves for several values of Nd 2 . The numerical errors are responsible for the noisy character 
of the simulation data for x > 10, however they are practically negligible for smaller values of x, the error bar being 
smaller than the symbol size. 
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X=t/Z 
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X=t/X 

FIG. 5. a) The results of the preceeding figures are reported in a plot of k(x). Only 200 points are selected in the range 
< x < 20. Crosses, plus signs and filled circles correspond to 7Vd 2 =0.064, 0.032 and 0.016, respectively. The continuous lines 
correspond to the fits obtained using Eq. (28) with the same proportionality factor, a — 0.424, between xo and Nd 2 in the 
three cases. The error bars are smaller than the symbol size, b) The region x > 5 is shown with a larger scale. At this new 
scale one can see the fluctuations of the numerical data giving an idea of the magnitude of the error bars (which are larger in 
the case Nd 2 = 0.016 due to the poorer statistics). 

Because of the high precision of our numerical simulation, we are able to verify the theoretical prediction concerning 
the role of the corridor effect. For this purpose we try to fit the numerical data to the function k(x) given by Eq. 
(B8), which does not take care of this effect, see Fig. 6. It is not possible to have a good fit with the same choice of 
the constant a as before. Thus we chose a different value, a = 0.283, to fit the maxima of the numerical curves. As 
it can be seen from Fig. 5, though the resulting fit could be considered as satisfactory, the difference is well beyond 
the numerical errors, in particular the numerical data points for 5 < x < 10 are systematically above the theoretical 
curves. Comparing Figs. 5 and 6, one can see that taking the corridor effect into account makes the agreement 
between the theoretical and numerical curves substantially better. 

Finally, in Fig. 7 we present the numerical results for the correction to the diffusion coefficient SD/D in units of 
Nd 2 /2tt versus \n{l/Nd 2 ), see Eq. (27). Data for both the hard-disk scattering and isotropic scattering are presented 
for several values of Nd 2 . For the case of isotropic scattering we get an excellent agreement with Eq. (27), in which 
the cutoff parameter is chosen as in Fig. 4: xo = 0A24Nd 2 . The slope of the dependence on \n(l/Nd 2 ) in the general 
case is given by Eq. (24). For the hard-disk scattering the slope is equal to 1.95 ± 0.10, which is close to what follows 
from Eq. (24). Indeed, for a given diameter d, the ratio of a tr a{Tx) for the two cases is equal to 2tt/3 = 2.094. 
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FIG. 6. a) The same numerical results are fitted by Eq. (B8), which does not take care of the corridor effect. Here, to fit the 
maximum values, we take a different proportionality factor than in figure 5, namely a — 0.283. b) The region x > 5 is shown 
with a larger scale. 
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FIG. 7. Numerical results for the diffusion coefficient as a plot of -(2n/Nd 2 )SD/D versus \n(l /Nd 2 ) for the hard disk 
scattering (open circles) and for isotropic scattering (filled circles). The linear behavior predicted by formula (27) in the 
isotropic case, with a = 0.424 as in Fig. 5, is represented by the straight line (of slope I). A linear fit through the three points 
in the hard disk case gives a slope of 1.95 ± 0.10 in good agreement with the theoretical prediction of 2.094. 



7. Conclusion 

In this paper, we have presented a simple approach to the problem of classical corrections to the Boltzmann 
equation, which are due to memory effects. We have performed both analytical and numerical calculations of the 
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velocity autocorrelation function for non-interacting particles scattered by randomly located centers in two dimensions. 

In the particular case of isotropic scattering, we were able to provide a full analytical expression for the correction 
to the Boltzmann result, due to returns of the particles to previously visited regions, in excellent agreement with the 
results of numerical simulations. Furthermore, the comparison between analytical and numerical results demonstrates 
that it is essential to take care of the "corridor effect" associated with backscattering events. In the case of hard 
disk scattering (Lorentz model), we obtain a good agreement for the dependence of the diffusion coefficient on the 
concentration of scattering centers. 

A full theory taking into account both classical and quantum (weak localization) memory effects is needed to 
understand the relative role of classical and quantum corrections depending on the ratio of the De Broglie wavelength 
to the scattering diameter. 

Appendix A 

For the case of scattering by hard disks, the Liouville equation has the form [6,7,15]: 

%+v%+v(f-+f + )f = 0, (Al) 
where the operators T_ and T + describe scattering by the disks. They are defined by the equations: 

f_/(r, <j>, t) = V d<j>'5((v - Ti - d/2)<r(0 - cf>')f(r, <f>, t), (A2) 



f+/(r, 0, t) = - V r #'<5((r - n + d/2)a(<j> - 0')/(r, ^ , t). (A3) 
i J o 

Here cr(0) = (d/4) sin((/>/2) is the differential cross-section for a disk of diameter d, the vector d is expressed through 
the velocities before and after the collision, v and v' by the relation (see Fig. 8): 

v' — v 

d = d ; (AA) 

^2(v 2 - vv') 

Equation (Al) is exact, the operators T_ and T + describing the mechanics of collisions (assuming that the disks do 
not overlap). The term T_/ describes the outflux of particles from the state with velocity v at the point + d/2 on 
the disk edge, while the term T + f describes the influx of particles to the state with velocity v at the point — d/2 
(see Fig. 8). 
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If we make the assumption that the distribution function f(r,<p,t) can be regarded as constant on distances on 
the order of d, we can neglect the vector d/2 in the arguments of the <5-functions in Eqs. (A2) and (A3). Then we 
immediately obtain Eq. (1). 



Appendix B 

For the case of isotropic scattering, we write the the Boltzmann equation for the function G in the form 



^ + ^ + 7(G - Go) = 6(r)5(<f>)5(t), (Bl) 



where 



G {T,t) = ±J G(r, 

is the isotropic part of G(r, <f), t). 

Performing the Fourier transformation of this equation in variables r and t, we obtain: 

- kv)G(k, 0, lu) + 7 (G(k, 0, lu) - G (k, lu)) = 5(<f>). (52) 

If we find G(k, <j>, lu) from this equation and take its average over the angle <f>, we will have an algebraic equation 
for Go(k, lu). Solving it and using Eq. (B2), we obtain the following expression for the function G(k,<f>,u): 

y-Jl m * + 

27r(7 — iu + ikvcos(<p — 4>k))(l — + ikv cos 4>k)(g(k, lu) — 7) ' 

where (f>k is the angle of the vector k and g(k, lu) = \/( 7 — iia) 2 + k 2 v 2 . 
Now we can easily find the expression for Gi(k, lu): 

Gl ( k)W ) = I 1 7exp(-^ fc ) g(k,Lu)- 1 + im 

' 27r(7 — iu> + ikv cos <j>k) 2mkv (7 — iu> + ikv cos <f>k){g(k,u>) — 7) 

The function Gi(0,t) entering Eqs. (18) and (19) is given by the Fourier transform of Eq. (B4) for r = 0. The 
first term in Eq. (B4) then becomes exp(— jt)5{v t). This contribution should be ommited, since it does not describe 
returns for t > 0. We begin the Fourier transformation of the second term by calculating the integral over the angle 
4>k using the formula [16] 

exp(—i<p) cos0 2n \f\ — a 2 — 1 

-d(p = 2 / -dtp — 

j_ n 1 + a cos q 11.. 

Then we have 



ff 1 + a cos <j> J 1 + a cos 4> a y/l — a 2 

Ol( o.o . r r * <#4^4. (*» 

(2tt)^ ./ x 2tt Jo 



g(k,uj){g(k,Lj) - 7) 

Instead of fc we now introduce a new variable z — g(k, lu) — 7 + iw, which gives us 

7 r dw _ lwt r 



GlM --JmwJ_ 00 to e ~*° Jo (z + 2( 7 -iu ) ))(z-iu 1 y {m) 
The integral over lu can be now easily calculated and the remaining integral over z gives: 

G i(M = Sh? (r^T^ + + 2cx P(-^)( Ei W - Ei(2x))), (Bl) 

where x = -ft. After replacing exp(— x) by exp(— x/2) in the first term to account for the corridor effect and introducing 
the cutoff at to — Z0/1 as described in Sect. 4, we obtain Eq. (26). 

If we neglect the corridor effect (but still introduce the cutoff at to), Eq. (B7) can be used directly to calculate the 
correction to the correlation function. Inserting Eq. (B7) into Eq. (18) we get, instead of k(x) given by Eq. (28), the 
following expression for k(x), which does not take into account the corridor effect: 

k(x) = -l + 2e- 2x (Ei(2x)-Ei{x)) +e- x 1 + (x - 2) (Ei(x) - C - In a;) + 2{x - 1) In 2 + (x + x Q ) hi(l + ^-) . (B8) 
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